Cosmological Constant: A Lesson from Bose-Einstein Condensates 
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The cosmological constant is one of the most pressing problems in modern physics. We address 
this issue from an emergent gravity standpoint, by using an analogue gravity model. Indeed, the 
dynamics of the emergent metric in a Bose-Einstein condensate can be described by a Poisson-like 
equation with a vacuum source term reminiscent of a cosmological constant. The direct computation 
of this term shows that in emergent gravity scenarios this constant may be naturally much smaller 
than the naive ground-state energy of the emergent effective field theory. This suggests that a 
proper computation of the cosmological constant would require a detailed understanding about how 
Einstein equations emerge from the full microscopic quantum theory. In this light, the cosmological 
constant appears as a decisive test bench for any quantum or emergent gravity scenario. 
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The cosmological constant [1] has been one of the most 
mysterious and fascinating objects for both cosmologist 
and theoretical physicists since its introduction almost 
a century ago [2]. Once called by Einstein his greatest 
blunder, it seems nowadays the driving force behind the 
current accelerated expansion of the universe. The expla- 
nation of its origin is considered one of the most funda- 
mental issues for our comprehension of general relativity 
(GR) and quantum field theory. 

Since this constant appears in Einstein equations as 
a source term present even in the absence of matter and 
with the symmetries of the vacuum (T^ v tx g^ v ), it is usu- 
ally interpreted as a "vacuum energy" . Unfortunately, 
this has originated the so-called "worst prediction" of 
physics. In fact, the estimated value, which is naively 
obtained by integrating the zero-point energies of modes 
of quantum fields below Planck energy, is about 120 or- 
ders of magnitude larger than the measured value. De- 
spite the large number of attempts (most notably super- 
symmetry [3, 4], which, however, must be broken at low 
energy) this problem is still open. We can summarize 
the situation by saying that, given the absence of cus- 
todial symmetries protecting the cosmological term from 
large renormalization effects, the only option we have to 
explain observations is fine tuning [5, 6]. 

This huge discrepancy is plausibly due to the use of 
effective field theory (EFT) calculations for a quantity 
which can be computed only within a full quantum the- 
ory of gravity (see, however, [7] for a proposal in the 
semiclassical gravity limit). Unfortunately, to date, we 
do not have any conclusive theory at our disposal. How- 
ever, the possibility of a failure of our EFT-based intu- 
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ition is supported by what can be learned from analogue 
models of gravity [8], given that, in these models, the way 
in which the structure of the spacetime emerges from the 
microscopic theory is fully under control. In [9, 10] it 
was shown that a naive computation of the ground-state 
energy using the EFT (the analogue that one would do 
to compute the cosmological constant), would produce a 
wrong result. The unique way to compute the correct 
value seems to use the full microscopic theory. 

Given the deep difference in the structure of the equa- 
tions of fluid dynamics and GR (and other gravitational 
theories), an accurate analogy cannot be performed at 
the dynamical level: indeed, this is forbidden by the ab- 
sence of diffeomorphism invariance and of local Lorentz 
invariance. However, in [11] it was shown for the first 
time that the evolution of part of the acoustic metric 
in a Bose-Einstein condensate (BEC) is described by a 
Poisson equation for a nonrclativistic gravitational field, 
thus realizing a (partial) dynamical analogy with Newto- 
nian gravity. Noticeably, this equation is endowed with a 
source term that is naturally identified as a cosmological 
constant, being there even in the absence of real phonons. 

In this Letter we will consider such analogue model for 
gravity and directly show that the cosmological constant 
term cannot be computed through the standard EFT ap- 
proach, confirming the conjecture of [9]. However, we 
find that also the total ground-state energy of the con- 
densate does not give the correct result: indeed, the cos- 
mological constant is comparable with that fraction of 
the ground-state energy corresponding to the quantum 
depletion of the condensate, i.e., to the fraction of atoms 
inevitably occupying excited states of the single particle 
Hamiltonian. In conclusion, the origin and value of such 
term teach us some interesting lessons about the cosmo- 
logical constant in emergent gravity scenarios. 

Settings. — The model used in [11] is a modified BEC, 
including a soft breaking of the U(l) symmetry associ- 
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ated with the conservation of particle number. This un- 
usual choice is a simple trick to give mass to quasipar- 
ticles that are otherwise massless by Goldstone's theo- 
rem. In second quantization, such a system is described 
by a canonical field \&T , satisfying [w(t, x), ^(t,* 1 )] — 
5 3 (x — x'), whose dynamics is generated by the grand- 
canonical Hamiltonian H = H — pN, where 
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and N is the standard number operator for In order 
for the interaction between bosons to be described by 
H, the gas must be dilute, i.e., pa 3 <C 1, where p is the 
density and a = Airgm/ti 2 is the s-wave scattering length. 
For more details on this model and on possible physical 
realizations, see [11, 12]. See also [13] for a generalization 
to condensates with many components. 

We describe the formation of a BEC at low tempera- 
ture through a complex function for the condensate 
and an operator cj> for the perturbations on top of it [14] : 



<3> = # (I- 
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Clearly, this is only an approximate characterization of 
the many body ground state. The validity of the mean 
field approximation must be checked, a posteriori, by 
controlling that the fluctuation (4> 2 ) is much smaller than 
l^o 1 2 = Po- If this is n °t so, the description of the effec- 
tive dynamics (e.g. the existence of an acoustic geometry 
where phonons propagate) does not hold any more. The 
canonical commutation relation for W implies 

-J-^x-x'). (3) 
Po(x) 
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We adopt the notation of [15], where a rigorous quanti- 
zation and mode analysis of the field <f> is presented for 
a standard BEC. Those results are here summarized and 
generalized to the [/(l)-breaking case of [11]. 

For a stationary condensate, d t ^o = and Eqs. (1) 
and (2) lead to a modified Gross-Pitaevski equation 
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For the aim of this Letter, it is enough to consider only 
homogeneous backgrounds. Thus, one can assume that 
V = and the condensate is at rest, such that \&o has 
a constant phase. For stability reasons, VPo must be real 
(^o* = = y/Po), and Eq. (4) simplifies to p = gpo — X. 

The equation for the quasiparticles is solved via Bo- 
goliubov transformation involving the Fourier expansion 
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where fik and at are quasiparticles' operators and the fac- 
tor */ po(27r) 3 has been inserted such that the Bogoliubov 
coefficients and t>k obey the standard normalization 
\ u k\ 2 — \vw\ 2 — 1- The dispersion relation is 
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describing massive phonons with ultraviolet corrections, 
mass M., and speed of sound c s [11] 
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As shown in [11], when the wavelength is larger than 
the healing length £ = h/mc Sl phonons propagate in 
an acoustic geometry with effective local Lorentz in- 
variance and their dispersion relation (6) is relativistic 
(quadratic). When k > the quartic term of the dis- 
persion relation (6) is instead dominant and the effective 
geometry is not defined. The Lorentz breaking scale Llv 
is thus identified with £. 

Standard manipulations give u^ 2 = (1 — D 2 .) ~ l and 
«k 2 = £*k M k) wrlere itk and i>k are chosen to be real and 
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Vacuum expectation values. — We can now compute the 
vacuum expectation value of H in the ground state |fl), 
the Fock vacuum of the quasiparticles (a k |f2) = 0, Vk). 
To this aim, it is convenient to expand H in powers of 
cf>: H ~ rlo + Hi + Hi, where Ho, Hi, and Hi con- 
tain, respectively, no power of cj>, only first powers, and 
only second powers, and higher order terms associated 
with quasiparticles' self-interactions are neglected. The 
energy density ho of the condensate (density of Ho) and 
the density hi of the expectation value of Hi are 
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while the expectation value of Hi vanishes because it 
contains only odd powers of Sk and a^. The integral in 
Eq. (9) is computed by using the above given expression 
for Vk. Applying standard regularization techniques [16] 
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where Fh is plotted in Fig. 1 (dashed line) and Fh(0)=l. 

The total grand-canonical energy density is therefore 
2 
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and it coincides with the well known Lee-Huang- Yang 
formula [17] when A = 0. 
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The number density operator N is analogously ex- 
panded in powers of 0: N — No + Ni + N 2 ■ The density 
of N Q is Po = |* | 2 , (JVi)n = 0, and p 2 = (N 2 )n is 



P2 = Pa 
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where satisfies F p (0) = 1 (see Fig. 1, dotted line). 
This is the number density of noncondensed atoms (de- 
pletion) and it is basically the magnitude of the fluctu- 
ations around the mean field. Note that poa 3 <^ lj a s 
described after Eq. (1). 

Furthermore, when A = 0, inverting the expression for 
total particle d ensit y, p = po + p 2 , one obtains, up to the 
first order in \ pa? 



Pa = P 
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which is the density of condensed atoms in terms of the 
total density p and the scattering length a [17]. In this 
case, p = gpo, such that the energy density e (density of 
(H) n = (n + pN) n ) is 
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This is the well known Lee-Huang- Yang [17] formula for 
the ground-state energy in a condensate at zero temper- 
ature. In general, when the U(l) breaking term is small, 
this term is expected to be the dominant contribution to 
the ground-state energy of the condensate. 

Analogue cosmological constant. — When the homoge- 
neous condensate background is perturbed by small in- 
homogeneities, the Hamiltonian for the quasiparticles can 
be written as (see [11]) 
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-ffquasip. is the nonrelativistic Hamiltonian for particles of 




FIG. 1: F h [dashed line, Eq. (10)], F p [dotted line, Eq. (12)], 
F M [dot-dashed line, Eq. (21)], and F A [solid line, Eq. (22)]. 



mass A4 [see Eq. (7)] in a gravitational potential 
(gpo + 3\)(gp + A) 
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and u(x) = [(po( x )/Poo) — l]/2, where p^ is the asymp- 
totic density of the condensate. Moreover, the dynamics 
of the potential <J> 3 is described by a Poisson-like equation 
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which is the equation for a nonrelativistic short-range 
field with length scale L and gravitational constant G^: 
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Despite the obvious difference between <j> 3 and the usual 
Newtonian gravitational potential, we insist in calling it 
the Newtonian potential because it enters the acoustic 
metric exactly as the Newtonian potential enters the met- 
ric tensor in the Newtonian limit of GR. The appearance 
of a short-range interaction in Eq.(17) is an artifact of 
the model. In [13] it has been shown how to obtain a 
long range analogue gravitational potential in a spinor 
BEC. However, the reasoning is identical in all the other 
relevant aspects, and the key result is unchanged. 

The source term in Eq. (17) contains both the contri- 
bution of real phonons (playing the role of matter) 
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where |C) is some state of real phonons, as well as a cos- 
mological constant like term (present even in the absence 
of phonons/matter) 

Ca = 2gpo(gpo + 3X)(g P o + X) Re + 1 ^ 

(20) 

Note that the source term in the correct weak field ap- 
proximation of Einstein equations is 4itGn(p + 3p/c 2 ). 
For standard nonrelativistic matter, p/c 2 is usually neg- 
ligible with respect to p. However, it cannot be neglected 
for the cosmological constant, since p\/c 2 = —pa- As a 
consequence Ca = — 2c 2 A, where A would be the GR 
cosmological constant. From Eq. (12) and evaluating 
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where F^^O) = 1 (see Fig. 1, dot-dashed line), we obtain 
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where F A = (2F p + 3i^)/5 (see Fig. 1, solid line). 
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Let us now compare the value of A either with the 
ground-state grand-canonical energy density h [Eq. (11)], 
which in [9] was suggested as the correct vacuum energy 
corresponding to the cosmological constant, or with the 
ground-state energy density e of Eq. (14). Evidently, A 
does not correspond to either of them: even when tak- 
ing into account the correct behavior at small scales, the 
vacuum energy computed with the phonon EFT does not 
lead to the correct value of the cosmological constant ap- 
pearing in Eq. (17). Noticeably, since A is proportional 
to \J poa, 3 , it can even be arbitrarily smaller both than 
h and than e, if the condensate is very dilute. Further- 
more, A is proportional only to the subdominant second 
order correction of h or e, which is strictly related to the 
depletion [see Eq. (12)]. 

Fundamental scales. — Several scales show up in this 
system, in addition to the naive Planck scale computed 
by combining h and the emergent constants Gn and c s : 

L * = « (—) ^ (A>a 3 r l7 V (23) 

V \gpoJ 

For instance, the Lorentz-violation scale Llv — £ 
(poa 3 )~ 1 / 2 a differs from Lp, suggesting that the breaking 
of the Lorentz symmetry might be expected at scale much 
longer than the Planck length (energy much smaller than 
the Planck energy), since the ratio Llv/-^p oc (p a 3 ) -1 / 4 
increases with the diluteness of the condensate. 

Note that Llv scales with p a 3 exactly as the range 
of the gravitational force [see Eq. (18)], signaling that 
this model is too simple to correctly grasp all the desired 
features. However, in more complicated systems [13], 
this pathology can be cured, in the presence of suitable 
symmetries, leading to long range potentials. 

It is instructive to compare the energy density corre- 
sponding to A to the Planck energy density: 

4ttGn fiGjj £ P \gPoJ 

(24) 

The energy density associated with the analogue cosmo- 
logical constant is much smaller than the values com- 
puted from zero-point-energy calculations with a cutoff 
at the Planck scale. Indeed, the ratio between these two 
quantities is controlled by the diluteness parameter poa 3 . 

Final remarks. — Taken at face value, this relatively 
simple model displays too many crucial differences with 
any realistic theory of gravity to provide conclusive evi- 
dences. However, it displays an alternative path to the 
cosmological constant, from the perspective of a micro- 
scopic model. The analogue cosmological constant that 
we have discussed cannot be computed as the total zero- 
point energy of the condensed matter system, even when 
taking into account the natural cutoff coming from the 
knowledge of the microphysics [9] . In fact the value of A 
is related only to the (subleading) part of the zero-point 



energy proportional to the quantum depletion of the con- 
densate. This holds also in a spinor BEC model, since 
the reasoning there is absolutely identical. The virtue of 
the single BEC model is to display the key physical re- 
sult without obscuring it with unnecessary mathematical 
complications, without loss of generality. Interestingly, 
this result finds some support from arguments within 
loop quantum gravity models [18], suggesting a BCS en- 
ergy gap as a (conceptually rather different) origin for 
the cosmological constant. 

The implications for gravity are twofold. First, there 
could be no a priori reason why the cosmological con- 
stant should be computed as the zero-point energy of the 
system. More properly, its computation must inevitably 
pass through the derivation of Einstein equations emerg- 
ing from the underlying microscopic system. Second, the 
energy scale of A can be several orders of magnitude 
smaller than all the other energy scales for the presence 
of a very small number, nonperturbative in origin, which 
cannot be computed within the framework of an EFT 
dealing only with the emergent degrees of freedom (i.e., 
semiclassical gravity) . 

The model discussed in this Letter shows all this ex- 
plicitly: the energy scale of A is here lowered by the di- 
luteness parameter of the condensate. Furthermore, our 
analysis strongly supports a picture where gravity is a 
collective phenomenon in a prcgcomctric theory. In fact, 
the cosmological constant puzzle is elegantly solved in 
those scenarios. From an emergent gravity approach, the 
low energy effective action (and its renormalization group 
flow) is computed within a framework that has nothing 
to do with quantum field theories in curved spacetime. 
Indeed, if we interpreted the cosmological constant as 
a coupling constant controlling some self-interaction of 
the gravitational held, rather than as a vacuum energy, 
it would immediately follow that the explanation of its 
value (and of its properties under renormalization) would 
naturally sit outside the domain of semiclassical gravity. 

For instance, in a group field theory scenario (a gen- 
eralization to higher dimensions of matrix models for 
two dimensional quantum gravity [19]), it is transparent 
that the origin of the gravitational coupling constants 
has nothing to do with ideas like "vacuum energy" or 
statements like "energy gravitates" , because energy it- 
self is an emergent concept. Rather, the value of A is 
determined by the microphysics, and, most importantly, 
by the procedure to approach the continuum semiclassi- 
cal limit. In this respect, it is conceivable that the very 
notion of cosmological constant as a form of energy in- 
trinsic to the vacuum is ultimately misleading. To date, 
little is known about the macroscopic regime of models 
like group field theories, even though some preliminary 
steps have been recently done [20]. Nonetheless, ana- 
logue models elucidate in simple ways what is expected 
to happen and can suggest how to further develop in- 
vestigations in quantum gravity models. In this respect, 



the reasoning of this Letter sheds a totally different light 
on the cosmological constant problem, turning it from a 
failure of effective field theory to a question about the 
emergence of the spacetime. 
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